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(57) Abstract 



The invention relates to a method for calculating a 
discrete orthogonal transformation according to the DIT 
(decimation in time) procedure in predetermined intermediate 
steps. Said method comprises the following steps: a) 
reading the data from a paged memory; b) carrying out 
the intermediate step predetermined by the transformation; 
c) saving the resulting data in an intermediate memory; 
and d) writing page by page the resulting data from the 
intermediate memory into the paged memory. The discrete 
orthogonal transformations can be FFT, IFFT, DCT, IDCT 
and structurally similar transformations. 



(57) Zusammenfassung 

Ein Verfahren zum Berechnen einer orthogonalen diskreten Transformation nach dem DIT-Verfahren in vorgegebenen Zwischen- 
schritten weist die folgenden Schritte auf: a) Lesen der Daten aus einem seitenweise organisierten Speicher; b) Durchfuhren des durch 
die Transformation vorgegebenen Zwischenschritts; c) Speichern der resultierenden Daten in einen Zwischenspeicher; und d) seitenweises 
Schreiben der resultie rende n Daten aus dem Zwischenspeicher in den seitenweise organisierten Speicher. Als diskrete orthogonale Trans- 
formationen kommen FFT, IFFT, DCT, IDCT und strukturell ahnliche Transformationen in Frage. 
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Beschreibung 

Verfahren und Vorrichtung zum Berechnen einer diskreten 
orthogonalen Transformation wie FFT Oder IFFT 

Die Erfindung betrifft ein Verfahren und eine Schaltung zur 
Berechnung der Fast-Fourier-Transformation sowie der dazu in- 
versen Transformation, im folgenden als FFT und IFFT bezeich- 
net, oder ahnlich strukturierter Transf ormationen, d.h. dis- 
krete orthogonale Transf ormationen . 

Die FFT und IFFT sind in der Nachrichtentechni k aufierst wich- 
tige Transf ormationen, da sie beispielsweise die Beschreibung 
eines Sachverhalts im Zeitbereich, in eine solche im Frequenz- 
bereich transf ormieren und umgekehrt. 

In der digitalen Signalverarbeitung wird haufig eine sog. N- 
Punkt diskrete Fourier-Transformation, im folgenden als DFT 
bezeichnet, berechnet, die wie folgt definiert ist: 

x < k ) = Z x(n > wnk k = 0, 1, , N-l 

mit W = e ("D ) 2n/N 

Die Komplexitat der Berechnung der DFT ist proportional zu 
0(N 2 ). Durch die Anwendung der FFT kann die Komplexitat der 
Berechnung auf 0{(Nlog(N)) reduziert werden. Dies geschieht 
durch eine hierarchische Aufteilung der Berechnung in Trans- 
f ormationen kurzerer Folgen. 

Zur Berechnung der FFT existieren zwei grundlegende Algorith- 
men. Der eine wird "Decimation in Frequency" (DIF) , der an- 
dere "Decimation in Time" (DIT ) genannt . Im folgenden wird 
beispielhaft der DIT-Algorithmus behandelt. 
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Zur Berechnung der FFT wird bevorzugt die sog. "in-place"-Va- 
riante verwendet, bei der errechnete Zwischenergebnisse der 
Butterf ly-Berechnung in denselben Speicher geschrieben wer- 
den, von wo sie wiederum gelesen und weiterverwendet werden, 
wie dies in der Fig. 1 dargestellt ist. Dadurch wird ein be- 
sonders sparsamer Umgang mit dem Speicher erreicht. 

In der Fig. 2 ist der Berechnungsvorgang der "in-place"-Vari- 
ante fur N = 8 in der Form eines Signal fluflgraphen darge- 
stellt. Wie aus der Fig. 2 hervorgeht, mussen die Daten am 
Anfang der Berechnung im Speicher in einer bestimmten Anord- 
nung vorliegen, die ublicherweise als "bit-reversed" bezeich- 
net wird. Am Ende der Berechnung kann das Ergebnis linear 
ausgelesen werden. Die Berechnung selbst geschieht in mehre- 
ren Stufen, wie es im Signalf luiigraphen der Fig. 2 darge- 
stellt ist. Im genannten Beispiel sind drei Stufen notwendig. 
Es werden jeweils 2 Daten aus dem Speicher ausgelesen, an- 
schlieiiend wird der Butterfly berechnet und die beiden Re- 
sultate werden wieder in dieselben Stellen im Speicher zu- 
rlickgeschrieben. Dabei liegen die Daten aber nicht unbedingt 
in benachbarten Speicherstellen . Ferner unterscheidet sich 
die Berechnung von einer Stufe zur nachsten. 

Wird die FFT in integrierten Schaltkreisen realisiert, so 
wird die Komplexitat vor allem durch den verwendeten Speicher 
bestimmt. Dabei sind groiie Speicher ublicherweise seitenweise 
aufgebaut, wodurch der Zugriff auf eine Speicherzelle inner- 
halb einer solchen Seite sehr schnell ist, beispielsweise ist 
ein derartiger Speicherzugrif f innerhalb eines Taktes durch- 
fuhrbar. Der Wechsel von einer Seite auf eine andere dauert 
aber wesentlich langer, d.h. mehrere Takte. Bei seitenori- 
entierten Speichern kann man den Durchsatz erhohen, indem man 
eine Seite eines Speichers moglichst vollstandig bearbeitet 
und erst dann wieder die Seite wechselt, weil Adressen der 
anderen Seite benotigt werden. Allerdings werden bei der oben 
genannten "in-place ,f -FFT die Daten prinzipiell sehr ungeord- 
net benotigt. Kleine Speicher haben diesen Nachteil nicht, da 
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hier der Zugriff auf die einzelnen Zellen uneingeschrankt 
moglich ist. 

Daher wird die Geschwindigkeit der FFT vor all em durch die 
Schnittstelle zum Speicher begrenzt. 

Der Erfindung liegt daher die Aufgabe zugrunde, ein Verfahren 
und eine Vorrichtung zur Berechnung diskreter orthogonaler 
Transf ormationen, insbesondere der FFT und I FFT, zu schaffen, 
womit eine schnellere Berechnung moglich ist. 

Die Aufgabe wird durch die Merkmale der unabhangigen Anspru- 
che gelost. Bevorzugte Ausgestaltungen der Erfindung sind Ge- 
genstand der Unteranspruche . 

Das erste er f indungsgemalie Verfahren zum Berechnen einer or- 
thogonalen diskreten Transformation nach dem DIT-Ver f ahren in 
vorgegebenen Zwischenschritten, wobei unter einem Zwischen- 
schritt die Addition sowie Multiplikation einer Verarbei- 
tungsstufe zu verstehen ist, umfafit die folgenden Schritte: 

a) Lesen der Daten aus einem seitenweise organisierten 
Speicher, 

b) Durchfuhren des durch die Transformation vorgegebenen 
Z wis chens chritts, 

c) Speichern der resul tierenden Daten in einem Zwischen- 
speicher, und 

d) Seitenweises Schreiben der resultierenden Daten aus dem 
Zwischenspeicher in den seitenweise organisierten Speicher. 

Das erf indungsgemalie Verfahren zum Berechnen einer orthogona- 
len diskreten Transformation nach dem DI F- Verfahren in vorge- 
gebenen Zwischenschritten weist die folgenden Schritte auf: 

a) Lesen der Daten aus einem seitenweise organisierten 
Speicher, 

b) Speichern der Daten in einem Zwischenspeicher, 

c) Durchfiihren des durch die Transformation vorgegebenen 
Zwischens chritts, und 
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d) Seitenweises Schreiben der resultierenden Daten aus dem 
Zwischenspeicher in den seitenweise organisierten Speicher. 

♦ 

Ein drittes er f indungsgemafies Verfahren zum Berechnen einer 
5 orthogonalen diskreten Transformation in vorgegebenen Zwi- ^ 
schenschritten weist die folgenden Schritte auf: 

a) Lesen der Daten aus zwei seitenweise organisierten Spei- 
cher, so daft das Lesen seitenweise organisiert ist, 

b) Durchfuhren des durch die Transformation vorgegebenen 
10 Zwischenschritts, und 

c) Seitenweises Schreiben der resultierenden Daten in die 
beiden seitenweise organisierten Speicher. 

Als diskrete orthogonale Transformation kommen beispielsweise 
15 eine FFT, I FFT , DCT oder IDCT sowie schematisch ahnlich orga- 
nisierte Trans formationen in Frage . 

Vorzugsweise weisen die Transformation eine identische Geome- 
trie fur j ede Stufe auf, was die Adressierung der Resultate 
20 erleichtert. Beispielsweise ist diese Bedingung fur eine FFT 
Oder I FFT nach Singleton erfullt. 

Entsprechend den obigen Verfahren weisen die er f indungsgema- 
Ben Vorrichtungen 

25 einen seitenweise organisierten Speicher, eine Rechen- 

einheit und einen nach der Recheneinheit angeordneten direkt 
organisierten Speicher, oder 

einen seitenweise organisierten Speicher, eine Rechen- 
einheit und einen vor der Recheneinheit angeordneten direkt 

30 organisierten Zwischenspeicher, oder 

zwei seitenweise organisierten Speicher und eine Rechenein- 
heit auf. . 

Dabei realisiert die Recheneinheit eine sog. Butterfly. ir 

35 



Bevorzugte Ausf uhrungsf ormen der Erfindung werden nachfolgend 
anhand der Zeichnungen erlautert. 
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Fig. 1 zeigt das Schema der Berechnung einer "in place"-FFT, 

Fig. 2 zeigt den Signalf luligraphen einer "in place"-FFT fur 
5 N = 8, 

Fig. 3 zeigt den Signalf luligraphen einer " Singleton"-FFT fur 
N = 8, 

10 Fig* 4 zeigt eine erste Ausf iihrungsf orm zur Berechnung einer 
FFT, und 

Fig. 5 zeigt eine zweite Ausf iihrungsf orm zur Berechnung einer 
FFT nach dem DIT-Verf ahren . 

15 

Fig. 1 zeigt die oben bereits erwahnte iibliche "in-place"-Be- 
rechnung einer FFT, wobei aus einem Speicher 1, ublicherweise 
eine DRAM ein Datenpaar gelesen wird, in einer Butterf ly-Ein- 
heit 2 verkniipft und wieder in den Speicher 1 geschrieben 
2 0 werden. 

Fig. 2 zeigt den grundlegenden DIT-Algorithmus einer FFT-Be- 
rechnung fur N = 8. Die Daten x(0) bis x(7) mussen zur Be- 
rechnung am Anfang "bit-reversed" im Speicher vorliegen. Dar- 

25 gestellt ist der Signalf lufigraph zur Berechnung der 8 Fou- 
rierkomponenten X(0) - X(7). Der Algorithmus sowie der Si- 
gnalf luiigraph ist A. V. Oppenheim, R.W.Schafer: "Digital Signal 
Processing", Prentice-Hall Inc., Englewood Cliffs, New Jer- 
sey, USA, 1975, S. 285 ff, beschrieben, so dafi hier nicht na- 

30 her darauf eingegangen werden braucht. Hinsichtlich der No- 
menklatur werden die Faktoren W°, W 1 und W 2 auch als Twiddle- 
Faktoren bezeichnet. 

Fig. 3 zeigt einen FFT-Algorithmus mit einer veranderten Ab- 
35 folge der Rechenoperationen, den sog. Singleton-Algorithmus , 
der in der. oben genannten Literaturstelle Oppenheimer-Schaf er 
auf S. 301 beschreiben ist, und der in dem erf indungsgemaJien 
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Verfahren verwendet wird. Es ist deutlich zu erkennen, daJJ 
die Daten sehr regelmafiig abgearbeitet werden, allerdings 
wird keine "in place"- Verarbeitung durchgef iihrt . Es wird da 
her der doppelte Speicherplatz im Vergleich zur "in place"- 
Berechnung benotigt . 

Vorteilhaf terweise ist die Berechnung in jeder Stufe vollig 
gleich aufgebaut, mit Ausnahme der Twiddle-Faktoren . D.h. di 
Stufen des beispielhaf ten Algorithmus einer FFT fur N = 8 
weisen eine identische Geometrie auf . Nach Abarbeitung einer 
Stufe wird vom zuletzt geschriebenen Speicher gelesen und in 
einen anderen Speicher geschrieben. 

Der Speicher kann linear gelesen werden. Ein lineares Schrei- 
ben ist bei dem Algorithmus jedoch noch nicht moglich, da im- 
mer ein Ergebnis in der oberen Halfte (voile Linie) und ein 
Ergebnis in der unteren Halfte ( gestrichelte Linie) gespei- 
chert werden muB . Allerdings kann in der oberen bzw. unteren 
Halfte linear geschrieben werden. 

Fig. 4 zeigt eine schematische Darstellung einer ersten er- 
f indungsgemaJJen Aus f uhrungsf orm zur schnellen Berechnung ei- 
nes Algorithmus, der ein ahnliches oder gleiches Schema wie 
der in Fig. 3 dargestellt Singleton-Algorithmus aufweist. 
Eine derartige Vorrichtung, die sowohl zur DIF- als auch DIT- 
Berechnung geeignet ist, ist der zur Berechnung notwendige 
Speicher in zwei gleich grofie seitenorientierte Speicher 3, 4 
unterteilt, die als DRAMs ausgelegt sind (DRAM = Dynamischer 
RAM) . Durch die Verwendung von zwei sei tenor ientierten Spei- 
chern 1 und 2 kann jeder Speicher in sich linear geschrieben 
werden, d.h. es wird jeweils ein Speicher fur die obere und 
untere Halfte des Algorithmus nach Fig. 3 verwendet so daii 
die Anzahl der langsamen Sei tenwechsel gering ist. Allerdings 
verdoppelt sich die Anzahl der Speicher und diejenige der 
Schnittstellen . 
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Eine, in der Fig. 5 dargestellt zweite Ausf uhrungsf orm zur 
schnellen Berechnung einer DIT-Transf ormation wie die der 
Fig. 3 ist die Verwendung eines kleinen schnellen Speichers 
5, der nach der Butterfly 2 angeordnet ist. Dort werden ei- 
nige Zwischenergebnisse der Berechnung zwischengespeichert , 
urn sie dann ohne standigen Seitenwechsel in den seitenorien- 
tierten Speicher 1 zu schreiben. Die Berechnung verlauft nun, 
indem aus dem Speicher 1 zwei Daten gelesen, der Butterfly 
zugefuhrt und die Zwischenergebnisse im schnellen SRAM-Spei- 
cher 5 (SRAM = Statisches RAM) gespeichert werden. Die Zwi- 
schenergebnisse werden dann seitenweise in den seitenorien- 
tierten Speicher 1 geschrieben. 

Bei einer Berechnung nach dem DIF-Verf ahren sitzt der stati- 
sche schnelle Speicher 5 am Eingang der Butterfly 2 (nicht 
dargestellt) . Ansonsten ist die Funktionsweise analog zu der- 
jenigen der erlauterten DIT-Variante . 

Die folgende Tabelle 1 soil das Adressierungsschema der Be- 
rechnung mittels Zwischenspeicherung verdeutlichen : 



Tabelle 1 



DRAM (RD) 


Butterfly 


SRAM (WR) 


SRAM (RD) 


DRAM (WR) 


0 


0 


0 






1 


16 


8 






2 


1 


1 






3 


17 


9 






4 


2 


2 


0 


0 


5 


18 


10 


1 


1 


6 


3 


3 


2 


2 


7 


19 


11 


3 


3 


8 


4 


4 


8 


16 


9 


20 


12 


9 


17 


10 


5 


5 


10 


18 


11 


21 


13 


11 


19 


12 


6 


6 


4 


4 
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13 


22 


14 


5 


5 


14 


7 


7 


6 


6 


15 


23 


15 


7 


7 


16 


8 


0 


12 


20 


17 


24 


8 


13 


21 


18 


9 


1 


14 


22 


19 


25 


9 


15 


23 


20 


10 


2 


0 


8 


21 


26 


10 


1 


9 


22 


11 


3 


2 


10 


23 


27 


11 


3 


11 


24 


12 


4 


8 


24 


25 


28 


12 


9 


25 


26 


13 


5 


10 


26 


27 


29 


13 


11 


27 


28 


14 


6 


4 


12 


29 


30 


14 


5 


13 


30 


15 


7 


6 


14 


31 


31 


15 


7 


15 








12 


28 








13 


29 








14 


30 








15 


31 



10 



In dem Beispiel hat der seitenorientierte Speicher eine Tiefe 
von N = 32, so dafi eine FFT mit 32 Punkten berechnet werden 
kann. Eine Seite des langsamen seitenorientierten Speichers 
habe beispielsweise P = 4 Adressen und der verwendete 
schnelle Speicher sei S = 4-P = 16 Adressen groB . Die Daten 
konnen linear in der Reihenfolge 0, 1, 2, ... gelesen werden. 
Beispielsweise wird zur Bearbeitung durch die Butterfly die 
Daten der DRAM Adressen 0 und 1 gelesen (RD) . Es ergeben sich 
durch die Zwischenergebnisse der Butterfly Adressen 0 und 16. 
Diese werden in die Adressen 0 und 8 des schnellen Zwischen- 
speichers .SRAM geschrieben (WR) . Mit einer Verzogerung, die 
zum teilweisen Fullen des Speichers benotigt wird, werden die 
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Daten seitenweise linear ausgelesen, also in der Tabelle 1 
unter SRAM (RD) die Adressen 0-3, anschliefiend 8-11, usw. Der 
Inhalt dieser SRAM-Adressen wird in die entsprechenden, durch 
die Transformation bedingten Adressen des langsamen sei- 
5 tenorientierten Speichers (Spalte DRAM (WR) ) seitenweise li- 
near geschrieben, also im Beispiel nach 0-3, 16-19, 4-7 usw. 
Durch die Verwendung des kleinen schnellen Speichers wird da- 
her der Zugriff auf den seitenorientierten Speicher DRAM der- 
art optimiert, dafi moglichst wenige langsarne Seitenwechsel 
10 notwendig sind. 

Die obengenannten und erlauterten Transf ormationen sowie 
Speichergroften dienen nur zur Erlauterung. In der Praxis ist 
der- langsarne _Speicher wesentlich grofier als der schnelle 

15 Speicher. Beispielsweise ist tiblicherweise der langsarne Spei- 
cher fur N = 8192 ausgelegt, wobei der langsarne Speicher eine 
Seitengroile von P = 16 aufweist. Daher hat der kleine Spei- 
cher eine GroBe von 64 Adressen. In einer integrierten Reali- 
sierung fallt der schnelle kleine Speicher daher f lachenmaJJig 

20 kaum ins Gewicht, allerdings wird die Berechnung der FFT oder 
IFFT oder ahnlicher diskreter orthogonaler Transf ormationen 
durch die Minimierung der Seitenwechsel des langsamen Spei- 
chers erheblich beschleunigt . 

2 5 Bezugszeichenliste 

1 seitenorientierter Speicher 

2 Butterfly 

3 seitenorientierter Speicher 
30 4 seitenorientierter Speicher 

5 schneller Speicher 
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Patentanspruche 

1. Verfahren zum Berechnen einer orthogonalen diskreten 
Transformation nach dem DIT-Verf ahren in vorgegebenen Zwi- 
5 schenschritten, 

dadurch gekennzeichnet, dali 

a) die Daten aus einem seitenweise organisierten Speicher 
(1) gelesen werden, 

b) der durch die Transformation vorgegebene Zwischenschritt 
0 durchgefiihrt wird, 

c) die resultierenden Daten in einen Zwischenspeicher (5) 
abgespeichert werden, und 

d) die resultierenden Daten aus dem Zwischenspeicher (5) 
seitenweise in den seitenweise organisierten Speicher (1) ge - 

5 schrieben werden. 

2. Verfahren zum Berechnen einer orthogonalen diskreten 
Transformation nach dem DIF-Verf ahren in vorgegebenen Zwi- 
schenschritten, 

dadurch gekennzeichnet, daft 

a) die Daten aus einem seitenweise organisierten Speicher 
(1) gelesen werden, 

b) die Daten in einen Zwischenspeicher (5) abgespeichert 
werden, 

O der durch die Transformation vorgegebene Zwischenschritt 
durchgefiihrt wird, und 

d) die resultierenden Daten seitenweise in den seitenweise 
organisierten Speicher (1) geschrieben werden. 

3. Verfahren zum Berechnen einer orthogonalen diskreten 
Transformation in vorgegebenen Zwischenschri tten, 

dadurch gekennzeichnet, da£ 
die Daten derart aus zwei seitenweise organisierten Speicher 
(3, 4) gelesen werden, daJ3 das Lesen seitenweise organisiert 
ist, 

der durch die Transformation vorgegebene Zwischenschritt 
durchgefiihrt wird, und 



) 
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die resultierenden Daten wieder seitenweise in die beiden 
seitenweise organisierten Speicher (3, 4) geschrieben werden. 

4. Verfahren nach einem der Anspruche 1 bis 3, dadurch 

5 gekennzeichnet, daft die diskrete orthogonale Transfor- 
mation durch eine FFT, IFFT, DCT Oder IDCT gebildet wird. 

5 . Verfahren nach Anspruch 4 , dadurch gekennzeich- 
net, daft die Transformation eine identische Geometrie fur 

10 jede Stufe aufweist. 

6 . Verfahren nach Anspruch 5, dadurch gekennzeich- 
net, daft eine FFT Oder IFFT nach Singleton verwendet wird. 

15 7. Vorrichtung zur Durchfiihrung des Verfahrens nach einem 
der Anspruche 1, 4-6, 

dadurch gekennzeichnet, daft 
die Vorrichtung einen seitenweise organisierten Speicher (1), 
eine Recheneinheit (2) und einen nach der Recheneinheit ange- 

20 ordneten direkt organisierten Speicher (5) aufweist. 

8. Vorrichtung zur Durchfiihrung des Verfahrens nach einem 
der Anspriiche 2, 4-6, 

dadurch gekennzeichnet, daft 
25 die Vorrichtung einen seitenweise organisierten Speicher (1), 
eine Recheneinheit (2) und einen vor der Recheneinheit ange- 
ordneten direkt organisierten Zwischenspeicher (5) aufweist. 

9. Vorrichtung nach einem der Anspruche 7 oder 8, da- 
30 durch gekennzeichnet, daft der seitenorientierte 

Speicher (1) im Verhaltnis zu dem direkt organisierten Zwi- 
schenspeicher (5) ein grofier Speicher ist. 

10. Vorrichtung nach Anspruch 9, dadurch gekenn- 

35 zeichnet, daft fiir den Zwischenspeicher (5) ein schneller 
Speicher verwendet wird. 
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11. Vorrichtung nach einem der Anspriiche 7- 10, d a d u r c h 
gekennzeichnet, daii der seitenorientierte Speicher (1) 
ein DRAM und der Zwischenspeicher (5) ein SRAM ist. 

12. Vorrichtung nach einem der Anspriiche 7 - 11, dadurch 
gekennzeichnet , daii der seitenorientierte Speicher ( 1 ) 
eine Grofie von 8 K Adressen und der Zwischenspeicher (5) eine 
Grofle von 32 - 64 Adressen hat. 



13. Vorrichtung zur Durchfuhrung des Verfahrens nach einem 
der Anspriiche 3, 4-6, 

dadurch gekennzeichnet, daii 
die Vorrichtung zwei seitenweise organisierten Speicher (3, 
4) und eine Recheneinheit (2) aufweist. 

14. Vorrichtung nach Anspruch 13, dadurch gekenn- 
zeichnet, daii die seitenorientierten Speicher (3, 4) 
gleich groll sind. 

15. Vorrichtung nach Anspruch 14, dadurch gekenn- 
zeichnet, daii der seitenorientierte Speicher (3, 4) eine 
Groiie von 4K Adressen aufweist. 



16. Vorrichtung nach 
gekennzeichnet, 
Butterfly realisiert . 



einem der Anspriiche 7 - 15, dadurch 
daii die Recheneinheit (2) durch einen 
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Fig. 3 
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